An Observational Method to Measure the Relative Fractions of Solenoidal and Compressible Modes in Interstellar Clouds

We introduce a new method for observationally estimating the fraction of momentum density (${\rho}{\mathbf{v}}$) power contained in solenoidal modes (for which $\nabla \cdot {\rho}{\mathbf{v}} = 0$) in molecular clouds. The method is successfully tested with numerical simulations of supersonic turbulence that produce the full range of possible solenoidal/compressible fractions. At present the method assumes statistical isotropy, and does not account for anisotropies caused by (e.g.) magnetic fields. We also introduce a framework for statistically describing density--velocity correlations in turbulent clouds.Comment: 20 pages, 13 figures, accepted for publication in MNRA

A Fluid-Dynamical Subgrid Scale Model for Highly Compressible Astrophysical Turbulence

We formulate and implement the Euler equations with SGS dynamics and provide numerical tests of an SGS turbulence energy model that predicts the turbulent pressure of unresolved velocity fluctuations and the rate of dissipation for highly compressible turbulence. We test closures for the turbulence energy cascade by filtering data from high-resolution simulations of forced isothermal and adiabatic turbulence. Optimal properties and an excellent correlation are found for a linear combination of the eddy-viscosity closure that is employed in LES of weakly compressible turbulence and a term that is non-linear in the Jacobian matrix of the velocity. Using this mixed closure, the SGS turbulence energy model is validated in LES of turbulence with stochastic forcing. It is found that the SGS model satisfies several important requirements: 1. The mean SGS turbulence energy follows a power law for varying grid scale. 2. The root mean square (RMS) Mach number of the unresolved velocity fluctuations is proportional to the RMS Mach number of the resolved turbulence, independent of the forcing. 3. The rate of dissipation and the turbulence energy flux are constant. Moreover, we discuss difficulties with direct estimates of the turbulent pressure and the dissipation rate on the basis of resolved flow quantities that have recently been proposed. In combination with the energy injection by stellar feedback and other unresolved processes, the proposed SGS model is applicable to a variety of problems in computational astrophysics. Computing the SGS turbulence energy, the treatment of star formation and stellar feedback in galaxy simulations can be improved. Further, we expect that the turbulent pressure on the grid scale affects the stability of gas against gravitational collapse.Comment: 19 pages, 16 figures, submitted to A&

An observational method to measure the relative fractions of solenoidal and compressible modes in interstellar clouds

This is the author accepted manuscript. The final version is available from Oxford University Press via the DOI in this record.We introduce a new method for observationally estimating the fraction of momentum density (ρv) power contained in solenoidal modes (for which ∇ · ρv = 0) in molecular clouds. The method is successfully tested with numerical simulations of supersonic turbulence that produce the full range of possible solenoidal/compressible fractions. At present, the method assumes statistical isotropy, and does not account for anisotropies caused by (e.g.) magnetic fields. We also introduce a framework for statistically describing density-velocity correlations in turbulent clouds. © 2014 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.A big thanks to Daniel Price for allowing us use of the auxiliary numerical simulations, to Maria Cunningham for allowing us access to the Delta Quadrant Survey data and to Dave Acreman for much-needed help with Figs 1 and 2. CB is funded in part by the UK Science and Technology Facilities Council grant ST/J001627/1 (‘From Molecular Clouds to Exoplanets’) and the ERC grant ERC-2011-StG_20101014 (‘LOCALSTAR’), both held at the University of Exeter. CF acknowledges funding provided by the Australian Research Council under the Discovery Projects scheme (grant DP110102191). Supercomputing time at the Leibniz Rechenzentrum (project pr32lo) and at the Forschungszentrum Jülich (project hhd20) are gratefully acknowledged. The software used in this work was in part developed by the DOE-supported ASC/Alliance Center for Astrophysical Thermonuclear Flashes at the University of Chicago

Implementing and comparing sink particles in AMR and SPH

We implemented sink particles in the Adaptive Mesh Refinement (AMR) code FLASH to model the gravitational collapse and accretion in turbulent molecular clouds and cores. Sink particles are frequently used to measure properties of star formation in numerical simulations, such as the star formation rate and efficiency, and the mass distribution of stars. We show that only using a density threshold for sink particle creation is insufficient in case of supersonic flows, because the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional physical collapse indicators have to be considered. We apply our AMR sink particle module to the formation of a star cluster, and compare it to a Smoothed Particle Hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas and sink particle properties between the AMR and SPH code.Comment: 6 pages, 4 figures, conference proceedings of IAU Symposium 270 (eds. Alves, Elmegreen, Girart, Trimble) simulation movies available at http://www.ita.uni-heidelberg.de/~chfeder/pubs/sinks/sinks.shtml?lang=e

Statistical analysis of the mass-to-flux ratio in turbulent cores: effects of magnetic field reversals and dynamo amplification

We study the mass-to-flux ratio (M/\Phi) of clumps and cores in simulations of supersonic, magnetohydrodynamical turbulence for different initial magnetic field strengths. We investigate whether the (M/\Phi)-ratio of core and envelope, R = (M/\Phi)_{core}/(M/\Phi)_{envelope} can be used to distinguish between theories of ambipolar diffusion and turbulence-regulated star formation. We analyse R for different Lines-of-Sight (LoS) in various sub-cubes of our simulation box. We find that, 1) the average and median values of |R| for different times and initial magnetic field strengths are typically greater, but close to unity, 2) the average and median values of |R| saturate at average values of |R| ~ 1 for smaller magnetic fields, 3) values of |R| < 1 for small magnetic fields in the envelope are caused by field reversals when turbulence twists the field lines such that field components in different directions average out. Finally, we propose two mechanisms for generating values |R| ~< 1 for the weak and strong magnetic field limit in the context of a turbulent model. First, in the weak field limit, the small-scale turbulent dynamo leads to a significantly increased flux in the core and we find |R| ~< 1. Second, in the strong field limit, field reversals in the envelope also lead to values |R| ~< 1. These reversals are less likely to occur in the core region where the velocity field is more coherent and the internal velocity dispersion is typically subsonic.Comment: 12 pages, 8 figures, accepted for publication in MNRA

Probing 3D Density and Velocity Fields of ISM in Centers of Galaxies with Future X-Ray Observations

Observations of bright and variable "reflected" X-ray emission from molecular clouds located within inner hundred parsec of our Galaxy have demonstrated that the central supermassive black hole, Sgr A*, experienced short and powerful flares in the past few hundred years. These flares offer a truly unique opportunity to determine 3D location of the illuminated clouds (with ~10 pc accuracy) and to reveal their internal structure (down to 0.1 pc scales). Short duration of the flare(s), combined with X-rays high penetration power and insensitivity of the reflection signal to thermo- and chemo-dynamical state of the gas, ensures that the provided diagnostics of the density and velocity fields is unbiased and almost free of the projection and opacity effects. Sharp and sensitive snapshots of molecular gas accessible with aid of future X-ray observatories featuring large collecting area and high angular (arcsec-level) and spectral (eV-level) resolution cryogenic bolometers will present invaluable information on properties of the supersonic turbulence inside the illuminated clouds, map their shear velocity field and allow cross-matching between X-ray data and velocity-resolved emission of various molecular species provided by ALMA and other ground-based facilities. This will highlight large and small-scale dynamics of the dense gas and help uncovering specifics of the ISM lifecycle and high-mass star formation under very extreme conditions of galactic centers. While the former is of particular importance for the SMBH feeding and triggering AGN feedback, the latter might be an excellent test case for star formation taking place in high-redshift galaxies.Comment: White paper submitted to the Astro2020 Decadal Surve

A method for reconstructing the variance of a 3D physical field from 2D observations: Application to turbulence in the ISM

We introduce and test an expression for calculating the variance of a physical field in three dimensions using only information contained in the two-dimensional projection of the field. The method is general but assumes statistical isotropy. To test the method we apply it to numerical simulations of hydrodynamic and magnetohydrodynamic turbulence in molecular clouds, and demonstrate that it can recover the 3D normalised density variance with ~10% accuracy if the assumption of isotropy is valid. We show that the assumption of isotropy breaks down at low sonic Mach number if the turbulence is sub-Alfvenic. Theoretical predictions suggest that the 3D density variance should increase proportionally to the square of the Mach number of the turbulence. Application of our method will allow this prediction to be tested observationally and therefore constrain a large body of analytic models of star formation that rely on it.Comment: 8 pages, 9 figures, accepted for publication in MNRA

Solenoidal versus compressive turbulence forcing

We analyze the statistics and star formation rate obtained in high-resolution numerical experiments of forced supersonic turbulence, and compare with observations. We concentrate on a systematic comparison of solenoidal (divergence-free) and compressive (curl-free) forcing, which are two limiting cases of turbulence driving. Our results show that for the same RMS Mach number, compressive forcing produces a three times larger standard deviation of the density probability distribution. When self-gravity is included in the models, the star formation rate is more than one order of magnitude higher for compressive forcing than for solenoidal forcing.Comment: 1 page, to appear in the proceedings of the IAU General Assembly Joint Discussion 14 "FIR2009: The ISM of Galaxies in the Far-Infrared and Sub-Millimetre", ed. M. Cunningha

Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH

We implemented sink particles in the adaptive mesh refinement (AMR) hydrodynamics code FLASH. Sink particles are created in regions of local gravitational collapse, and their trajectories and accretion can be followed over many dynamical times. We perform a series of tests including the time integration of circular and elliptical orbits, the collapse of a Bonnor-Ebert sphere and a rotating, fragmenting cloud core. We compare the collapse of a highly unstable singular isothermal sphere to the theory by Shu (1977), and show that the sink particle accretion rate is in excellent agreement with the theoretical prediction. To model eccentric orbits and close encounters of sink particles accurately, we show that a very small timestep is often required, for which we implemented subcycling of the N-body system. We emphasize that a sole density threshold for sink particle creation is insufficient in supersonic flows, if the density threshold is below the opacity limit. In that case, the density can exceed the threshold in strong shocks that do not necessarily lead to local collapse. Additional checks for bound state, gravitational potential minimum, Jeans instability and converging flows are absolutely necessary for a meaningful creation of sink particles. We apply our new sink particle module for FLASH to the formation of a stellar cluster, and compare to a smoothed particle hydrodynamics (SPH) code with sink particles. Our comparison shows encouraging agreement of gas properties, indicated by column density distributions and radial profiles, and of sink particle formation times and positions. We find excellent agreement in the number of sink particles formed, and in their accretion and mass distributions.Comment: 30 pages, 17 figures, ApJ accepted, simulation movies available at http://www.ita.uni-heidelberg.de/~chfeder/videos.shtml?lang=e